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ABSTRACT 

The computational study of water jet cooled hot side heat sink for thermoelectric refrigerator has been ana- 
lyzed to study the effects of change in nozzle length (l) and diameter (d) on various parameters like the convective heat 
transfer coefficient, thermoelectric resistance, pressure drop and Nusselt number. It was observed from the analysis that 
nozzle configuration ofl=3 mm and d=l mm gave maximum value ofthe heat transfer coefficients as compared to other 
nozzle geometries considered at various flow rates. Computational studies have been done to analyzed the heat sink per- 
formance by introducing dimple surfaces of three different diameters of lmm, 2mm and 3mm at the inlet port of 8 mm. 
Results were further more analyzed for nozzle configuration of 1=3 mm, d=lmm with respect to the temperature 
differences oftarget surface and exiting fluid, for various flow rates from 1-5 litres per minute. Results showed that heat 
sink model with dimpled surfaces of diameter 3 mm gives maximum Nusselt number and convective heat transfer coeffi- 
cient; high Pressure drop at lesser thermal resistance. 
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INTRODUCTION 

Thermoelectric refrigeration is the latest growing technology based on Peltier element that is environment 
friendly and in which heat transfer can take place through solid and liquid medium [1]. Its development tends to 
eliminate the problem of leakage of refrigerant that is generally found for the case of vapour compression refrige- 
ration system (VCRS) and also problem of bulky size like in case of vapour absorption refrigeration system 
(VARS). In Peltier element, there are a number of thermocouples arranged electrically in series and thermally in 
parallel through which when the electric current is passed, a temperature differential is generated across the two 
sides on the surface. The main advantage of using Peltier element is that it can be even installed in very small 
space of one cubic inch [1]. Thermoelectric refrigeration has a wide range of applications such as for medicine and 
electronics where precise temperature control is required, portable ice boxes, transport of perishable products, mil- 
itary, aerospace and many more. Thermoelectric refrigeration (TER) is not widely used due to its low coefficient 
of performance (COP around 0.38 - 0.45) in comparison to VARS (COP around 0.6 - 0.7) and VCRS (COP 
around 2.6 - 3.0) concluded by Riffat et al [2]. This is the reason it is not being considered for many of the prac- 
tical applications. The COP of TERs can be enhanced if the efficient heat sink is used. Until now, many types of 
heat sinks have been developed like air cooled heat sink, free and submerged jet liquid cooled heat sink and micro 
jet heat sink. Manoj et al. [3] reviewed the development in TER and TEAC systems and found that these systems 
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can be used as a substitute to traditional refrigeration and air-conditioning systems, but due to less development of these 
systems, the efficiency of these systems is around 5 to 15 % in comparison to conventional systems, which have efficiency 
around 40 to 60 %. Wei et al [4] investigated that Capability of TEAC and TEHS driven by PV/T modules and found out 
that the COP of TEAC is more than 0.45. The COP of TERS and TEAC can be enhanced if the difference to the tempera- 
ture between hot side and cold side of Peltier element is low. Cosnier et al. [5] presented a numerical model of TEAC and 
TEHS and also investigated it experimentally and found that COP between 1.5 and 2.0 can be achieved by maintaining the 
temperature difference of 5 °C and by applying an electric current of 4A. Shen et al [6] analyzed air-conditioning system 
based on thermoelectric radiant with an electric supply of 1.2 A and obtained COP of 1.77. The system can also be used as 
an air heating system just by reversing the input current. Zhou and Yu [7] found that in designing of TERS, an important 
role is played by the hot side heat sink and also the thermal conductance of the heat sink for hot side should be optimized. 
Hamed et al. [8] found that the cold side performance of the thermoelectric modules can be enhanced by adjusting the tem- 
perature of hot side by water flow at a constant rate, and then it can also be used as an air cooling system. Riffat et al. [9] 
found efficiency of thermoelectric refrigerators can be enhanced with an optimum design of the heat sink for hot side. It 
has also been found that the COP of TERS in which water is used to cool hot side heat sink doubles than TERS with the air 
cooled hot side heat sink. The water cooled heat sink can be of two types, one is the free jet impingement heat sink and 
other is the submerged jet impingement heat sink. Robinson and Schnitzler [10] found in experimental investigation that 
the submerged jet impingement heat sink gives the higher heat transfer coefficient than that of the heat sink with free jet 
impingement for a given pumping power. Robinson [11] concluded that due to higher pressure drop in case of the micro jet 
heat sink, the submerged jet impingement heat sink is better as for practical considerations. Single submerged jet im- 
pingement has been compared with multiple submerged jet impingement heat sink by many researchers. Tie et al. [12] 
stated that later gives more uniformity of temperature distribution on target surface and enhancement in heat transfer coef- 
ficient. Experimentally and numerically, it has been found that many factors like jet to target distance, impingement veloci- 
ty, jet diameter, pressure drop, liquid flow rate effect the performance of the hot side heat sink. 


Table 1 


Nomenclature 


s w 

Jet to Wall Distance (mm) 

VCRS 

Vapour compression refrigeration system 

d 

Nozzle diameter (mm) 

VARS 

Vapour absorption refrigeration system 

1 

Nozzle Length (mm) 

TER 

Thermoelectric refrigeration system 

A t 

Exposed target area of heat source 
(mm 2 ) 

TEHS 

Thermoelectric heating system 

Q 

Heat transfer (W) 

TEAC 

Thermoelectric air conditioning 

H 

Convective heat transfer coeffi- 
cient (W/m 2 K) 

COP 

Coefficient of Performance 

AT 

Temperature difference between 
target surface and exiting fluid 
(K) 

H 

Jet to target distance (mm) 

Rth 

Thermoelectric Resistance (K/W) 

s 

Jet to jet distance (mm) 

lpm 

Litre per minute 


In the present paper, convective heat transfer coefficient, Nusselt number, temperature variations were analyzed 
for various flow ranges. Modifications were made in the hot side heat sink based on Karwa et al. [13]. 1 and d has been 
varied to three different values of lmm, 2mm and 3mm to analyze the effect on Convective heat transfer coefficient, ther- 
moelectric resistance (R t h), Nusselt number and Pressure drop. Computational analysis were performed for the model with 
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introduction of dimpled surfaces at the inlet to the pipe for nozzle configuration of l=3mm, d=lmm and comparative study 
were also made between the two model with/without dimpled surfaces. 

COMPUTATIONAL MODEL 

Multiple submerged jet impingement heat sink model was selected to study the performance of TERS. Deionized 
Water was used as cooling fluid. As shown in Figure l(a), the water enters from inlet port and then reaches the plenum 
(conical section). At the base of plenum there are numbers of nozzles, arranged in 5 rows x 5 columns. Jet to jet centre dis- 
tance (s) was given as 6.25mm and jet to outside wall distance (s w ) as 5.5mm. The walls in model were given thickness of 
2mm. It was observed that chamfering and contouring at inlet of nozzle decrease the pressure drop in case of submerged 
heat sink. A fillet of 0.2 mm radius was given at the inlet of nozzles. Then water passes through nozzles and strikes on hot 
side of peltier element (which is the target surface). After striking the target surface with high velocity, the water exits per- 
pendicular to the direction of target surface and then leaves from the model by outlet port. Interfacial thermal resistance 
was neglected as fluid directly strikes on target surface. There are two models, model (a) and (b), as shown in Figure 1. 
The model (a) is without dimples on which effect of change in nozzle length and nozzle diameter has been studied on heat 
transfer coefficient and other parameters. The model (b) is with dimpled surfaces. The effect of dimpled surfaces was stu- 
died on the temperature difference between target surface and exiting fluid by varying the inlet velocity. The various di- 
mensions of the heat sink (in mm) are shown in Figure l(b) and Figure 3. Length of the inlet/outlet port is 20mm, height of 
plenum is 15mm and the length of the connecting part from inlet port to plenum is also 20mm. The base of the heat sink 
model is 70x7Omm 2 . 



(a) Without Dimpled Surfaces (b) With Dimpled Surfaces 

Figure 1: Sectional View of the Presented Heat Sink Model 
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Water flow 



* Dimensions not to scale 

Figure 2: Direction of the Flow on the Peltier Element 


Table 2: Mesh Independency Test 


Solid Works 
Mesh 

Resolution 

No. of 

Total cell 
Elements 

No. of 

Solid cell 
Elements 

No. of 

Fluid cell 
Elements 

No. of Partial 
Cell Elements 

Convective Heat Transfer Coefficient, h (W/m 2 K) 






1 lpm 

2 lpm 

3 lpm 

4 lpm 

5 lpm 

3 

29022 

5915 

6800 

16307 

2979.2 

3291.12 

3445.32 

3538.82 

3583.3 

4 

70818 

19270 

21278 

30270 

2122.7 

2340.99 

2467.16 

2486.53 

2512.7 

5 

192238 

60790 

69964 

61484 

5455.79 

6330.26 

6772.25 

7040.76 

7230.6 

6 

192238 

60790 

69964 

61484 

5460.4 

6329.73 

6778.54 

7034.24 

7228.4 

7 

192238 

60790 

69964 

61484 

5465.17 

6332.22 

6772.16 

7033 

7229 


Figure 2 shows the position of peltier element (Target surface) along with the direction of the cooling fluid (wa- 
ter) and also various dimensions of the nozzle. Figure 3 shows the wire frame model of the selected heat sink and its di- 
mensions. 


20 



Figure 3: Dimensions of the Selected Heat Sink 


NUMERICAL PROCEDURE 


For the computational heat transfer and fluid flow analysis, Solid Works 2012 edition was used. A mesh inde- 
pendency study was performed as shown in Table 1. Total number of 192238 cells at mesh resolution of 5 was selected for 
computational analysis 
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Heat transfer coefficient was calculated by Fourier equation 

Q = hA t AT (1) 

Where, Q denotes heat transfer rate, h denotes average convective heat transfer coefficient, AT denotes the differ- 
ence of temperature between average temperature of target surface and bulk average temperature of the fluid, A t denotes 
target surface area which is 38x38 mm 2 . 

Nusselt number was calculated as 

Nu = (2) 

Where, Nu denotes the Nusselt number, k denotes thermal conductivity of flowing fluid (water). 1 is the characte- 
ristic length which was assumed as the diameter of the nozzle. 

Thermal resistance was calculated as 

R* = f (3) 

Where, R th denotes thermal resistance, AT denotes temperature difference between target surface and bulk average 
temperature of the fluid and Q denotes heat transfer rate. 

BOUNDARY EQUATIONS 

Volumetric flow rate was varied from 1 lpm to 5 lpm; accordingly, mass flow rate value was given as input. Out- 
let boundary condition was selected as environmental pressure. Constant heat flux rate of 34626 W/m 2 was given at target 
surface (50 W for 38 x 38 mm 2 ) where the water coming from the jet strikes. Inlet temperature of fluid was set as 293.20 K. 
Model walls were assumed to be adiabatic. Conduction and radiation effects were neglected. 

MODEL VALIDATION 

Present model has been validated against the pressure drop readings over a range of Re (Reynolds number) from 
847 to 4246 for nozzle of d = 1 mm. It can be observed in Figure 4 that present computational model has the same trend of 
readings as observed in similar heat sink model proposed by Karwa et al [13]. The present work experienced lower pres- 
sure drop (7.6%) and some variation in trend as compared to Karwa et al. [13] due to introduction of fillet of radius 0.2 mm 
at the inlet of nozzle. The boundary conditions and nozzle size in both simulations are selected same. 
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Figure 4: Model Validation by Comparing the Results of Present Computational 
Model with Karwa Et Al. [13] for Pressure Drop V/S Re 

RESULTS AND DISCUSSIONS 

Convective Heat Transfer Coefficient 

The Figure 5(a) shows the variation in convective heat transfer coefficient (h) with the flow rate. It increases with 
increase in fiow rate. The reason behind the increment of h can be considered due to mixing of striking fluid at fast rate. 
This decreases the temperature of target surface, but as the volumetric flow rate increases, the bulk average temperature of 
exiting fluid also decreases. In case of nozzle geometry l=3mm, d=lmm, which is the case of maximum value of convec- 
tive heat transfer coefficient (in the model without dimple surface). In this case, as the flow increases from 1 lpm to 5 lpm, 
with 1 lpm increase in each time, value of h increases by 16.84 %, 6.62 %, 3.77 % and 3.16% respectively. Maximum in- 
crement in h with flow rate was observed when flow rate value changed from 1 lpm to 2 lpm 



Flow Rate (Ipm) 


(a) for Various Nozzle Configurations without 
Dimpled Surfaces 



Figure 5: H V/S Flow Rate 

It was also observed from the above Figure 5(a) that h increases with increment in length of nozzle (1) and with 
decrement in nozzle diameter (d) for a fixed flow rate. Increment of h w.r.t. length variation is lesser as compare to incre- 
ment w.r.t. decrease in diameter. Maximum nozzle length assumed is 3mm and minimum nozzle diameter is 1 mm. there- 
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fore as we discussed, according to that maximum h value is observed in case of nozzle geometry l=3mm and d=lmm. min- 
imum value of h was observed at l=lmm and d=3mm. the value of h increases by 40.77 % when nozzle size was changed 
from 1=1 mm, d=3mm to l=3mm, d= 1 mm at llpm. if decrease the diameter from lmm to any lower value, there are chances 
of clogging of nozzles practically and nozzle length more than 3mm can also be restricted due to space constraint. 

As shown in Figure 5(b), by introducing dimpled surfaces in the present model, at flow rate of 1 lpm h value de- 
creases than that of without dimpled surface model, when dimple surface of lmm, 2mm and 3mm introduced respectively. 
At flow rate of 2 lpm with dimpled surfaces of diameter lmm, h value also decreases, but as dimpled surfaces of diameter 
2mm and 3mm used h value increased at the same flow rate. When it was increased more than 2 lpm h value increases 
when dimple surface of lmm, 2mm and 3mm used. Maximum increase in h was observed at flow rate of 4 lpm when dim- 
pled surfaces of 3mm were used. 

THERMAL RESISTANCE 

As shown in Figure 6(a) it was observed that the thermal resistance (R th ) decreases with the increase in flow rate. 
Decrement in R th follows the trend of polynomial of 2nd order (with the regression coefficient around 0.95-0.98). The rea- 
son for the decrease of R th with the increase of flow rate is attributed due to more fluid component striking the target sur- 
face thus leading to less temperature difference AT (using Rth = AT/Q). At nozzle geometry l=lmm, d=3mm when flow 
rate was varied from 1 lpm to 5 lpm, the decrement of Rth was observed to be 21.60- 3.58 %. 



(a) for Different Nozzle Configurations Without 
Dimpled Surfaces 
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(b) for l=3mm,d=lmm, With Dimpled Surfaces 


Figure 6: Thermal Resistance V/S Flow Rate 

It was observed that with increase in flow rate from 1 lpm to 51pm, at constant length of lmm, when diameter of 
nozzle increase from 1 mm to 3mm, the thermal resistance increases by 14.32-33.33% and at constant diameter of 1 mm, 
when length of the nozzle increase from lmm to 3mm, it decreases by 3.32-5.29%. At constant heat flux and length of 
nozzle as the diameter increases R th increases, this may be attributed due to less striking velocity of the fluid on the target 
surface thereby decreasing the convective coefficient which leads to increase in AT. The increase in R th with an increase in 
d is much more than that of decrease in 1. Minimum thermal resistance is observed at nozzle geometry l=3mm, d=lmm at 
flow rate of 5 lpm and maximum thermal resistance at l=lmm, d=3mm with flow rate of 1 lpm. R th decreases by 15.43% 
when nozzle size is changed from 1=1 mm, d=3mm to l=3mm, d=lmm at flow rate of 51pm. 

As shown in Figure 6(b) when dimpled surfaces were introduced in the heat sink model, R th was increased. With 
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increase in diameter of dimpled surfaces up to 3mm, R th decreases, but not much decrement was observed. Diameter of 
dimpled surfaces can be increased more than 3mm,to observe the variation in R th . 

PRESSURE DROP 

It was observed from the Figure 7(a) that as flow rate increases, the Pressure drop (dP) also increase and it follows 
the polynomial trend (of 2 nd order). At a given flow rate and fixed nozzle length, dP was observed to increase with decrease 
in diameter of the nozzle. This increment in dP with decrease in d at a given flow rate can be attributed due to lesser the 
diameter at the nozzle, more will be the obstruction in the flowing fluid. It was also observed to increase with increase in 
nozzle length at fixed flow rate and nozzle diameter. The increment in dP v/s flow rate with an increase of 1 is observed to 
be at a higher rate when d=lmm. Maximum pressure drop is observed at nozzle geometry of l=3mm, d=lmm with flow 
rate of 5 lpm. Minimum pressure drop was observed at nozzle geometry of l=lmm, d=3mm with flow rate of 1 lpm. At this 
nozzle geometry, minimum pressure observed was 0.29 kPa when flow rate was llpm and maximum pressure drop ob- 
served was 6.73 kPa when flow rate was 5 lpm. When nozzle geometry was changed from 1=1 mm, d=3mm to l=3mm, 
d=lmm at 51pm, dP increased to 46.92 kPa.At constant flow rate, nozzle length and diameter (l=3mm and d=lmm), with 
the introduction of dimpled surfaces, dp was observed to be increased by 0.5-2.5% when dimples diameter changed from 1 
mm to 3mm. Maximum pressure drop of 47.37kPa was observed with dimple surface of 3mm at flow rate of 5 lpm. 
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Figure 7: Dp Y/S Flow Rate 


NUSSELT NUMBER 


Nusselt number (Nu) tells us about the enhancement in heat transfer rate when we give motion to stagnate fluid. 
As shown in Figure 8(a) Nu was observed to increase with the increase in flow rate. Its variation followed the logarithmic 
trend (with the regression coefficients around 0.99). It was also observed that Nu was increased with the increase in 1 at 
constant flow rate and d. Increase in Nu with increase in d was also observed at constant flow rate and 1. At a constant 
l=3mm and flow rate of llpm, when the nozzle diameter was changed from l-3mm, the percentage increase in Nu was ob- 
served as 129%. At a constant nozzle diameter of lmm and flow rate of llpm when the nozzle length was changed from 1- 
3mm, the percentage increase in Nu was observed as 5.58%. Thereby it was analyzed that Nu increases more rapidly with 
the nozzle diameter than that with nozzle length. Maximum Nu was observed at nozzle geometry of l=3mm, d=3mm. 
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As shown in Figure 8(b), by introducing the dimpled surfaces in the heat sink model, not much variation was ob- 
served in Nu. When flow rate was more than 21pm, it was observed that Nu increases with increase in diameter of dimpled 
surfaces. Maximum value of Nu was found for dimpled surfaces of diameter 3mm. 




(a) for Different Nozzle Configurations Without (b) for l=3mm,d=lmm, With Dimpled Surfaces 
Dimpled Surfaces 

Figure 8: Nu V/S Flow Rate with Change in Nozzle Geometry (Without Dimpled Surfaces) 

TEMPERATURE DIFFERENCE BETWEEN TARGET SURFACE AND EXITING FLUID 



1 2 3 4 5 

Flow Rate (Ipm) 


Figure 9: dT ( Target Surface - Fluid) v/s Flow Rate with Change in Nozzle Geometry 


Temperature difference between target surface and exiting fluid depends on various parameters, including but not 
limited to flow rate of cooling liquid, inlet fluid temperature and nozzle geometry. It was analyzed in the model without 
dimpled surfaces that nozzle geometry of l=3mm and d=lmm gave lower thermal resistance, and higher heat transfer coef- 
ficient. At this nozzle geometry, the pressure drop was also high, which can be reduced if the lesser flow rate may be used. 
Dimpled surfaces were introduced at the inlet section of the model for this nozzle geometry.The results of four configura- 
tions of the model were observed which are shown in Figure 9. The first configuration was the heat sink model without 
dimpled surface; in that situation, the temperature difference between the target surface and exiting fluid was lowest at 
flow rate of 1 lpm. The second configuration was the model with dimpled surfaces of diameter lmm, in which the tempera- 
ture difference was observed to be highest at same flow rate of 1 lpm. The third configuration was the model with dimpled 
surfaces of diameter 2mm, in which the temperature difference observed was closer to the model with dimpled surface of 
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diameter d=3mm (which is the fourth configuration of the model) at flow rate of llpm. With the further increment of flow 
rate to 2 lpm and more, the temperature difference was observed to be lesser in case of the model with dimpled surfaces of 
3mm diameter. Minimum temperature difference was observed for this geometry at flow rate of 4 lpm. 

CONCLUSIONS 

The simulation study of water jet cooled submerged jet hot side heat sink has been done to analyze the effect of 
change in nozzle length and diameter on various parameters (Convective heat transfer coefficient, Thermal resistance, 
pressure drop, Nusselt number). With the study of these parameters, it has been observed that 

• Maximum convective heat transfer coefficient was observed for the model with dimpled surfaces of diameter 
3mm at flow rate of 5 lpm. 

• Minimum thermal resistance was observed in the same conditions for maximum heat transfer coefficient. 

• Minimum pressure drop was found in the model without dimpled surfaces at nozzle geometry of 1=1 mm, d=3mm 
with flow rate of 1 lpm. 

• Maximum Nusselt number was observed for model without dimpled surfaces at nozzle geometry of l=3mm, 
d=3mm with flow rate of 5 lpm. 

• Minimum temperature difference between target surface and exiting fluid was observed for the model with dim- 
pled surfaces of diameter 3mm at flow rate of 4 lpm. 

In this heat sink model, jet pattern (square array) can be modified to get the more uniform temperature on target 
surface, dimpled surfaced of higher diameter (than 3mm) can also be designed to study the effect on heat sink model para- 
meters and turbulence effect may be studied. Thermoelectric heat sink model’ s plenum design can be modified to enhance 
the model’ s performance. 
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